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FORWARD ERROR CORRECTOR 
CROSS-REFERENCE TO RELATED APPLICATION 

This patent application is a continuation of U.S. Patent 
Application Serial No. 10/382,400 filed March 6, 2003, which is 
5 a continuation of U.S. Patent Application Serial No. 09/951,998 
filed September 12, 2001 now U.S. Patent No. 6,539,516 issued 
March 25, 2003, which is a continuation of U.S. Patent 
Application Serial No. 09/437,448 filed November 9, 1999 now 
U.S. Patent No. 6,317,858 issued November 13, 2001, which claims 
10 priority on the basis of the provisional patent application 
Serial No. 60/107,879 filed November 9, 1998. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

15 The present invention relates to an apparatus for 

correcting errors present in stored or transmitted data; and, 
more particularly, to an apparatus for evaluating an error 
evaluator polynomial, an error locator polynomial and a 
differential polynomial which are used in correcting errors in 

2 0 the data encoded by using an algebraic code, such as a 
Reed-Solomon code. 

2 . Description of Related Art 

Noise occurring during a process of transmitting, storing 
2 5 or retrieving data can in turn cause errors in the transmitted, 
stored or retrieved data. Accordingly, various encoding 
techniques, having the capability of rectifying such errors, for 
encoding the data to be transmitted or stored have been 
developed. 

30 In such encoding techniques, a set of check bits is 

appended to a group of message or information bits to form a 
codeword. The check bits, which are determined by an encoder, 
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are used to detect and correct the errors. In this regard, the 
encoder essentially treats the bits comprising the message bits 
as coefficients of a binary message polynomial and derives the 
check bits by multiplying the message polynomial R(x) with a 
5 code generator polynomial G(x) or dividing R(x) by G(x) , to 
thereby provide a codeword polynomial C(x). The code generator 
polynomial is selected to impart desired properties to a 
codeword upon which it operates so that the codeword will belong 
to a particular class of error-correcting binary group codes 

10 (see, e.g., S. Lin et al . , "Error Control Coding: Fundamentals 
and Applications", Prentice-Hall, 1983). 

One class of error correcting codes is the well-known BCH 
(Bose-Chaudhuri-Hocquenghen) codes, which include the 
Reed-Solomon ("RS") code. The mathematical basis of the RS code 

15 is explained in, e.g., the aforementioned reference by Lin et 
al. and also in Berlekamp, "Algebraic Coding Theory" , 
McGraw-Hill, 1968, which is further referred to in U.S. Pat. No. 
4,162,480 issued to Berlekamp. The aforementioned references 
are hereby incorporated by reference in pertinent part. 

20 

SUMMARY OF THE INVENTION 

The invention herein provides a method and apparatus for 
decoding an algebraic-coded message. The method can include the 
steps of determining a discrepancy indicator, with the 

25 discrepancy being between a calculated and a predicted value; 
determining an error locator polynomial using a selected class 
of error correction algorithms, such as, for example, a 
Berlekamp-Massey algorithm; and detecting an uncorrectable 
message using the selected error correction algorithm. The 

30 apparatus is composed of storage devices which can include 
recirculating storage devices; arithmetic components attached to 
the storage devices, the components operating over a Galois 
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Field on selected contents of the storage devices; and an 
uncorrectable message detector, connected with the storage 
devices and the arithmetic components. 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is an illustration of a algebraic decoder 
according to the invention herein; 

Figure 2 is a data flow diagram of a modified Berlekamp- 
Massey algorithm according to the invention herein; 
10 Figure 3 is a block diagram illustrative of an exemplary 

embodiment of the present invention; 

Figure 4 is a block diagram of a circular syndrome 
generator according to the present invention; and 

Figure 5 is a block logic diagram of a logic register which 
15 can be used in the circular syndrome generator illustrated in 
Figure 4 . 

EXEMPLARY EMBODIMENTS OF THE INVENTION 

The invention herein provides an apparatus for and a method 

2 0 of decoding algebraic codes, including BCH codes, and more 
specifically, Reed-Solomon codes, such that uncorrectable 
messages, or portions of received encoded data, are detected. 
Furthermore, the invention herein provides for a more area- 
efficient device implementation of the aforementioned method. 

25 For the purposes of illustration, the present invention will be 
described in terms of a subset of the BCH codes, namely Reed- 
Solomon (RS) codes. 

The Reed Solomon (RS) encoding technique appends to each 
block of k user data symbols 2t redundancy symbols to create an 

30 encoded message block (where t represents the designed symbol 
error correcting capacity of the code) . These 2t symbols, or 
elements, are selected from the Galois Field to be the roots of 
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the generator polynomial. Therefore, there are k+2t symbols in a 
RS-encoded message block. The entire message block is viewed as 
a polynomial and evaluated as a polynomial at some Galois Field 
element. The Galois Field element at which the polynomial is 
5 evaluated will be located at one roots of the generator 
polynomial that are used to create the RS code. The RS code 
views the n-bit symbols as elements of a Galois Field (GF ( 2n) ) . 
A Galois field is a finite field, the elements of which may be 
represented as polynomials in a, where a is a root of an 

10 irreducible polynomial of degree n. The RS codeword consists of 
a block of n-bit symbols. Typically, n = 8 and the 8-bit symbols 
are referred to as bytes. Constructing the Galois field GF(2n) 
requires a defining polynomial F(x) of degree n. In addition, a 
primitive element p is chosen so that every nonzero element of 

15 GF(2n) is a power of (3. The element (3 is not necessarily a root 
of F(x) . 

A RS codeword C is viewed as a polynomial C(x) and the 
redundancy symbols are chosen so that the roots of C(x) include 
the roots of a generator polynomial G(x) whose roots are 2t 

20 consecutive powers of (3. The k user data symbols are viewed as 
the high order coefficients of a degree k+2t-l polynomial, and 
the redundancy symbols are the coefficients of the remainder 
when this polynomial is divided by G(x) . 

The process of corrupting the original code block C(x) with 

25 errors can be viewed as adding an error polynomial E(x) to C(x). 
The resultant corrupted polynomial is known as the received 
polynomial R(x) , where R (x) =C (x) +E (x) . The v non-zero terms of 
the error polynomial contain the necessary information to 
completely reconstruct the original data C(x), since each term 

30 corresponds to a symbol error location and magnitude. 

Typically, RS decoding is a tripartite analysis: (1) 
syndrome computation; (2) solution of the error magnitude and 
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locator polynomials; and (3) error location and magnitude 
estimation by respective implementations of, for example, a 
Chien search and the Forney algorithm. The syndromes contain 
error information divorced form the actual information that is 
5 intended to be analyzed for errors. The error locator 

polynomial provides information regarding the location of an 
error in the received signal, and the magnitude of the error can 
be determined by using both the magnitude and the locator 
polynomials. 

10 The thrust of the RS error correction procedure is to 

reconstruct the error polynomial E(x). Three polynomials are 
used to correct a received polynomial R(x) : S (x) , a syndrome 
polynomial; A(x), an error locator (or error location) 
polynomial; and M(x) an error magnitude polynomial. The 

15 syndromes are computed by evaluating the polynomial R(x) at all 
roots of G(x) . These values are called syndromes and the 
syndrome polynomial S(x) has these values as coefficients. The 
syndrome polynomial S (x) is used to determine the existence of 
errors. The error locator polynomial A(x) and the error 

20 magnitude polynomial M(x) are computed from S(x) by a key 
equation solver. The roots of the error locator polynomial A(x) 
indicate positions in the data that are erroneous and both the 
error locator polynomial A(x) and the error magnitude polynomial 
M(x) are used to determine the true values of the erroneous 

2 5 data. 

Two frequently-used RS error correction algorithms are the 
Berlekamp-Massey and the Euclid algorithms. The present 
invention recasts the Berlekamp-Massey algorithm such that the 
inversion process typically associated with the traditional 
30 Berlekamp-Massey (tBM) algorithm is eliminated. This is 

important because the inversion process includes determining the 
reciprocal of certain Galois field elements using division. 
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Division is a time consuming arithmetic operation, the 
implementation of which can occupy needed component area in a 
device design. Therefore, the present invention can be 

particularly advantageous where area-efficient layout of a 
5 decoder device is desirable. 

For further elaboration of the decoding process over Galois 
fields, including tBM, Chien searching, and Forney's Algorithm, 
see Theory and Practice of Error Control Codes by Richard E. 
Blahut (Addison-Wesley, 1983) which is incorporated by reference 

10 in pertinent part herein. 

Figure 1 illustrates an implementation of this algorithm, 
in which a raw received signal 1 is directed to RS decoder unit 
2 that is used to determine the error locations and error 
values. Signal 1 is provided to syndrome generator 3 and delay 

15 unit 4. In syndrome generator 3, the several syndromes 5 
associated with the selected encoding are derived and 
transmitted to polynomial solver 6. The syndrome generator 3 
calculates one syndrome for each of the 2t roots of G(x). 
Polynomial solver 6 utilizes the syndromes to determine the 

2 0 coefficients of the error location polynomial A(x) 7 and the 
coefficients of the error magnitude polynomial M(x) 8, which in 
turn are transmitted to error estimator 9. Estimator 9 
calculates error signal 10 which is combined in summer 11 with 
delayed raw received input 12 to provide corrected data 13. 

25 Estimator 9 can include Chien search unit 14 which utilizes the 
error location polynomial A(x) to search for the roots of the 
error locator polynomial, r lf . . . , r v . Typically, the Chien 
search unit 14 employs a root finding technique which involves 
evaluating the error locator polynomial at all elements in the 

30 field GF(2 n ). The roots of the error locator polynomial r lf . . . 
, r v determine the error locations. The error values are then 
determined using Forney's algorithm unit 15. The delayed raw 



-6- 



51318/RJP/B600 

received input 12 is then corrected using the output of the 
Forney algorithm unit 15 and the raw received input which is 
transmitted by delay unit 4. 

Traditionally, the Berlekamp-Massey (tBM) algorithm, which 
5 usually is realized in polynomial solver 6 can described by: 

A r = Z A ( i r-1) Sr-j (1) 

j=o J 

Lr = 5r(r-Lr-l + (l-6r)L r -l (2) 



10 
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r=l,..., 2t where 6 r = 1 if both A r * 0 and 2L r _i < r-1, and 
otherwise 5r = 0 . Then A (2t) (x)is the smallest-degree polynomial 
with the properties that A 0 (2t) =l, and 

15 

S r + n ZAj 2t) Sr-j = 0 r = L 2 t + 1— 2t 

j=l J 

where initial conditions are A (0) (x)=l, B (0) (x)=l ; and L 0 =0. 
It is evident that the inversion indicated in Eq. 3 
requires a division operation. 

The tBM algorithm is capable of properly decoding messages 
that can be decoded properly, however if there is an 
uncorrectable case which is detectable as being uncorrectable, 
the uncorrectable error may be undetected and the message 
decoded as if it did contain a correctable error. Many times, 
this improperly decoded message can create additional 
difficulties because the error may propagate through other 
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processes in the system which employs tBM. 

According to the invention herein, the modified Berlekamp- 
Massey (mBM) can be described by the following equations: 



A i =lAt 1 1 S i " j (4) 



Ai = A_Ai-i + xAiBi-i (5) 



Bi = 



'M-i 
xB^ 



A_=A; 



(6a) 



10 (6b) 
where: A * 0 
B 0 = 1 
A_!= 1 

15 Utilization of mBM for RS decoding can be advantageous 

because: (1) inversion is eliminated; (2) the control structure 
associated with the mBM algorithm is simplified relative to that 
of tBM; and (3) the termination conditions of tBM are modified 
such that if the code is uncorrectable, errors otherwise 

20 undetected by tBM, are detected and flagged as such. 

One implementation of the mBM algorithm is as follows, as 
represented in Pascal code: 



PROCEDURE FindLocatorBMC ( VAR Syndrome , Locator : Polynomial ; VAR OK: BOOLEAN ); 

25 

VAR Cnt: 0 . .MaxParitySyms-1 ; { Loop Index } 

Pwr: 0 . .MaxParitySyms ; { Power Counter } 

State: (Alpha, Beta) ; { State Machine State } 

Deg: INTEGER; { Degree } 

30 

Del, DelO: Words; { Discrepancies } 

J: INTEGER; { Del Index } 
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45 



TempPoly : Polynomial ; 
B : Polynomial; 

BEGIN 

B.L. := 0; B.D [0] :=1; 
Locator.L := 0; Locator. D [0] :=1, 
Deg : = 0; Pwr : = 0; DelO := 1; 
State := Alpha; 



{ Temporary Polynomial} 
{ Secondary Polynomial} 

{ BEGIN FindLocatorBMC} 
{ Initial B } 
{ Initial Locator Poly } 
{ Cntr Initialization } 
{ Machine State } 



FOR Cnt := Pari tySyms-1 DOWNTO 0 DO BEGIN { Algorithm Loop 

Del := 0; { Calculate Del 

FOR J := 0 TO LOCATOR.L DO 
IF Syndrome. L >= ( ParitySyms-l-Cnt-J) THEN 

Del:= Add( Del, Multiply{ Locator . D [J] , Syndrome . D [ Pari tySyms-1 -Cnt- J] )) ; 



TempPoly : = 



{ Do Common Update 



} 



PolyAdd( WordTimes( Locator, DelO ), PolyShift( WordTimes ( B, Del ) , 1) ) 



IF (State=Alpha) AND (DeloO) THEN BEGIN 
{ writeln( stderr, ' B<-L' );} 
B := Locator; 
DelO := Del 
END 

ELSE BEGIN 
{ writeln( stderr, ' B<-xB' );} 
B := PolyShift( B, 1 ) 
END; 

IF State=Alpha THEN BEGIN 

IF Del=0 THEN Pwr := Pwr +1 

ELSE State := Beta 

END 

ELSE BEGIN 

Deg := Deg+1; 

IF Pwr = 0 THEN State := Alpha 

ELSE Pwr := Pwr-1 

END; 

Locator := TempPoly 
END; 

Locator := PolyDenormalize ( Locator, Deg) ; 

OK := State=Alpha 

END; 



{ Do Step A 

{ Do Step A 

{ Do Step B 

{ Do Step B 



} 
} 

} 

{ State is Alpha } 
{ Increment Power Cntr} 
{ Update Next State } 
{ State is Alpha } 



{ State is Beta 



} 



{ Update Next State } 
{ Decrement Power Cntr} 
{ State is Beta } 
{ Update Locator } 
{ Algorithm Loop } 

{Update Locator Degree} 

{ END FindLocatorBMC } 



50 



Often, when a forward error corrector properly detects an 
uncorrectable error, the existence of such an error usually is 
verified in a process by which the syndrome polynomials are 
recomputed. This approach can carry a substantial penalty 
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relative to the process efficiency. Instead, an embodiment of 
the invention herein, having an improved control structure, 
verifies the existence of an uncorrectable error by checking the 
state of polynomial solver 6 at the end of the polynomial 
5 solving process. 

Figure 2 exemplifies an embodiment of the process 2 0 
implementing the aforementioned improved control structure in 
the context of the mBM algorithm recited in Equations 4, 5, and 
6 (a) -(b). Although the implementations described herein are 

10 postured for standard RS codes having a block length of, for 
example, 2 55 elements, such implementations also may be used in 
the context of extended Reed-Solomon codes which, in the example 
herein, would have 256 elements in the message block, i.e., have 
256 elements in associated the Galois Field. It is desirable 

15 that, in step 21, the control variables DEG, PWR, and STATE, as 
well as error locator variables be initialized. It further is 
desirable to iterate through steps 23, 24, 25, and 26, 2t times, 
where 2t is the number of syndrome polynomials to be evaluated, 
and t is the error correcting capability of the preselected 

20 code. Thus, at step 30, a counter tracking the number of 
completed iterations is employed. No additions or subtractions 
are needed in implementing the control variables, and only count 
up or down functions are used. Step 23 essentially implements 
Equation 4, in which the discrepancy value DEL, associated with 

25 a particular iteration, is determined. Similarly, step 24 
implements Equation 5 in which the error locator polynomial is 
updated. In step 25, auxiliary polynomial Bi is updated 
according to Equation 6a in substep 27, or Equation 6b in 
substep 28, based on conditions determined by logic 26. For 

30 logic 29, it is desirable for both STATE = ALPHA AND DEL < > 
zero to direct the data flow via an implementation of Equation 
6a in. substep 27; otherwise substep 28 is used, implementing 
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Equation 6b. Unlike the tBM algorithm where the polynomial 
shift term (1 — 5 r )x in Equation 3 has been normalized, the mBM 
algorithm does not require normalization, avoiding an 
inversion/division operation. After the auxiliary polynomial is 
5 updated in step 25, the controller state is updated in step 26. 

In general, the degree of the error locator polynomial is 
tracked by DEG, which is an upcounter descriptive of the true 
degree of the error locator polynomial and, thus, the number of 
errors in the message block. It also is desirable to construct 

10 an error locator polynomial who roots equate to the locations of 
an error. Essentially, process 20 attempts to synthesize a 
linear feedback shift register (LFSR) that predicts the values 
of the syndrome polynomial. Such a LFSR can be useful to find 
the error locations. Discrepancy value, DEL, then exposes a 

15 discrepancy between the predicted value of the syndrome 
polynomial, and the value of the current syndrome polynomial, 
and invites further processing to discover the location of the 
errors. PWR is a counter that keeps track of the number of 
times that the controller previously remained in STATE = ALPHA. 

2 0 It is desirable to have the STATE remain in control state BETA 

for a count equivalent to the number of times that STATE 
previously remained in control state ALPHA. 

For the purposes of the invention herein, STATE can be used 
to (1) determine whether the error correction analysis ought to 
25 follow the flow of Equation 6a or 6b; (2) assist in determining 
whether the degree of the error locator polynomial ought to be 
increased; and (3) whether there exists an uncorrectable error. 
At the end of 2t iterations, the value of STATE is once again 
determined at step 35. If the result is STATE = ALPHA, then the 

3 0 code is potentially valid; on the other hand, if STATE = BETA, 

then the error is flagged as uncorrectable. Potentially valid 
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codes where STATE = ALPHA at step 35, also can be subjected to 
additional validation before being subsequently decoded. 
Indeed, in one subsequent operation, the number of the error 
locator polynomial zeroes is compared with the value of DEG. A 
5 discrepancy between these two values also is indicative of an 
uncorrectable error. 

Figure 3 is an exemplary embodiment of a polynomial solver 
using the mBM algorithm. Solver 5 0 can include syndrome 
generator register 51, recirculating syndrome register 52, 
10 first delay register 53, locator polynomial (Aj)register 54, 
auxiliary polynomial (Bi) register 55, second delay register 56, 
first multiplier 57, second multiplier 58, adder 59, Aj register 

60, and A register 61. In another embodiment, syndrome 
generator 51 can be separate from solver 50. It is desirable 

15 for multipliers 57, 58, and adder 59 to operate on Galois Field 
elements. It also is desirable for register 51 to be logically 
arranged as a circular register or loop, such that particular 
register values can be used in a pre-defined sequence. 
Furthermore, it is desirable that registers 52, 54, and 55 be 

20 logically arranged as push-down stacks or FIFOs, and also that 
the values contained therein rotate synchronously. The error 

locator polynomial Aj are arranged in register 54 such that the 

least significant coefficient is at the top and the stack "grows 

down" as subsequent values are determined. It is desirable for 
2 5 the syndrome recirculating register 52 to operate such that the 

least significant coefficient is aligned with the bottom of the 

register, and the most significant with the top. 

In the example of Figure 3, the error correcting 

capability, t, is selected to be 5 elements and, thus, syndrome 
30 register 51 is designed to employ 2t, or 10, individual 

elements. Additionally, register 52 is selected to use t 
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elements, register 54 is chosen to employ t+1 elements, and 
register 55 is intended to operate with t elements. 

Initially, recirculating syndrome register 52 is pre-loaded 
with zeroes. As expressed in the aforementioned algorithm, an 
5 initial step involves calculating the discrepancy value DEL, 
which can be stored in register 60. A previous value for DEL, 
DELO, is provided in register 61. To calculate the initial value 
for DEL, first syndrome value S 0 is shifted into register 52, 
which value is received in first multiplier 57 along with the 

10 initial value of DEL, namely, DELO, and the t-th value in locator 
polynomial register 54. After the indicated multiplication and 
creation of a DEL value, the values in registers 52, 54, and 55 
are shifted down by one element. At first, .the values in 
register 52 are zero, however, with subsequent clocking, 

15 successive values of Si enter register 52 and are recirculated 
therethrough, for the clock cycles equivalent to i = 0 to 2t-l. 
As So exits register 52 into first multiplier 57, corresponding 

values of Aq are also transmitted to first multiplier 57 such 

that the value SqAq is determined. 

20 Concurrently with this calculation, value B 0 from register 

55 is multiplied with then extant value for DEL in second 

multiplier 58 and the result is summed with SqAq in adder 59 to 

produce the next value for DEL. This value of DEL is used to 
produce the next value for the error locator polynomial, namely, 

25 Aj . After this calculation, value So is recirculated such that 

it bypasses first delay register 53 and re-enters recirculating 
register 52 at the top of the stack during the next clock cycle. 
During this next clock cycle, syndrome value Si is aligned, and 

multiplied, with A 0 , and S 0 is aligned, and multiplied, with Aj . 

3 0 This process repeats such that each of the syndrome values 
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properly iterates through in the evaluation of the locator 
polynomial . 

With the above information, a skilled artisan would be able 
to see the manner in which the values for error locator 

5 polynomial Aj and auxiliary polynomial Bi are determined. Where 
it is desired to rotate the values of Bi through register 55, 
second delay register 56 is bypassed. On the other hand, where 
it is desirable to utilize a previously calculated value of Bi, 
as indicated in the aforementioned algorithm, the value of Bi is 
10 directed into second delay register 56. 

Continuing in Figure 3, calculation of the magnitude 
polynomial will be described. At the completion of 2t iterations 
as described . above, register 52 will contain values S 4 -S 8 . In 
essence, determination of the magnitude polynomial can be 

15 modeled as M(x) = A(x)S(x) mod X , in which the multiplication 

will be truncated after the 2t term. Indeed, only t terms need 
be determined under the assumption that no uncorrectable error 
was encountered. During the final iterations of the calculation 
of the error locator polynomial, register 52 is loaded with 
20 zeros such that, at the completion of 2t iterations, all storage 
locations in register 52 contain a zero value. After 2t 
iterations, the values in register 51 will be restored to their 
original positions, register 52 will contain all zero values and 

register 54 will contain the error locator polynomial, Aj . In 

2 5 am manner similar to the computation of the locator polynomial 
coefficients, the error magnitude coefficients are calculated 
iteratively. After t iterations, S 0 will be at the logical 
bottom of register 52 and Aq at the logical bottom of register 
54. At the completion of t+1 iterations, the product of 

30 multiplier 57 will be SqAq , the first term of the error 
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magnitude polynomial. The output of adder 59 is directed to the 
logical top of register 55, which now will be used to build the 
magnitude polynomial. After iteration t+2, syndrome value So 
will be aligned with locator value A± , giving the product SqAj ; 

5 syndrome value SI will be aligned with Aq , giving the product 
SiAo ' the summation of which giving S0A1 + S1A0' which is the 
second term of the error magnitude polynomial. This process will 
continue until all values of M(x) are so determined. At 
iteration 2t, all of the coefficients for the error magnitude 
10 polynomial will have been calculated. At this point, data flow 
of the error locator polynomial in register 54 and the error 
magnitude polynomial in register 55 can be directed out of 
solver 50. 

Figure 4 illustrates one embodiment of a circular syndrome 
15 generator 70 that can be employed as register 51 in Figure 3, 
modified to accommodate an error correcting capability of t=8. 
Figure 5 is an embodiment 72 of one of the individual registers 
71 in circular syndrome generator 70 in Figure 4. Although a 
circular syndrome generator is shown, it is by no means the only 
20 form of syndrome generator that can be employed as register 51. 

The foregoing merely illustrates the principles of the 
invention, and it will thus be appreciated that those skilled in 
the art will be able to devise various alternative arrangements 
which, although not explicitly described herein, embody the 
25 principles of the invention within the spirit and scope of the 
following claims. 
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